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Abstract 

Employing quantum kinetic equations we study the formation of binary correlations 
in plasma at short time scales. It is shown that this formation is much faster than 
dissipation due to collisions, in hot (dense) plasma the correlations form on the 
timescale of inverse plasma frequency (Fermi energy) . This hierarchy of characteris- 
tic times is used to derive analytical formulae for time dependency of the potential 
energy of binary interactions which measures the extent of correlations. We discuss 
the dynamical formation of screening and compare with the static screened result. 
Comparisons are made with molecular dynamic simulations. In the low tempera- 
ture limit we find an analytical expression for the formation of correlation which is 
general for any binary interaction. It can be applied in nuclear situations as well as 
dense metals. 



Recent lasers allow one to create a high density plasma within few femto sec- 
onds and observe its time evolution on a comparable scale [1,2]. Naturally, this 
plasma is highly excited at the beginning and relaxes towards equilibrium by 
various mechanisms that might be dominant at some stage and sub-dominant 
in another one. The best known regimes are the fast local equilibration of 
electron and hole distributions due to binary collisions, and the slow global 
relaxation via diffusion, recombination, dissipation of energy into the host 
crystal, etc. There are, however, even faster processes than the local equilibra- 
tion. These processes dominate during the very first stage of the relaxation, 
the so called transient regime. In this Letter we discuss the transient regime 
in terms of the energy balance. 
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The transient regime has been already discussed from many different angles 
and for various systems (classical - quantum, non-degenerate - degenerate). 
Let us briefly review some of these approaches. We will show that they are 
equivalent, at least with respect to the energy balance. 

Conceptually the simplest is the molecular dynamics. One takes N particles, 
distributes them randomly into a box and let them classically move under 
Coulomb forces due to their own charges. Those particles which are very close 
will be expelled from each other. Their first movement thus forms correlations 
which lower the Coulomb energy Vq = e 2 jr. This build up of screening stops 
when the effective Debye potential Vd — e 2 e~ Kr /r is reached. An important 
question is whether the long-range or the short-range charge fluctuations dom- 
inate in this process. In the former case, the characteristic time of the transient 
period should be the inverse plasma frequency r c = l/cu p . In the latter case 
we do not know. 

If one "measures" all states of the systems by their distance from equilibrium, 
instead of formation of correlations one has to talk about their decay. Our 
presumption that the transient period is appreciably shorter than the local 
relaxation is thus just Bogolyubov's principle of decay of correlations. The 
decay of correlations is linked with an alternative approach to the transient 
period, the formation of quasi-particles which has been numerically studied 
within Green's functions [3]. We note that very similar transient behavior 
has been observed for the nuclear matter [4,5], i.e., the formation/decay of 
correlation is a rather general phenomenon. 

The first concept is based on the two-particle space correlations while the 
second one on the single-particle excitations. The quantity which allows us 
to follow both pictures in a unified manner is the energy of the system. It 
is composed from the kinetic energy (j^j and the correlation energy -E cor r — 
|(^d) — (Yc)i where (Vc) subtracts the background. In lowest order interaction 
and classical limit the correlation energy takes the Debye Hiickel form 

E co „= 1 -lim[V D (r)-V c (r)] = -^f. (1) 

Z r— >0 Z 



Of course, the total energy conserves, 

/•-;•„,, /£- > -( — ). (2) 




where (Jm) * s ^ ne i n itia.l value of the kinetic energy. We will monitor the time 
dependency of the transfer of the correlation energy into the kinetic one. 

It is more convenient to calculate the kinetic energy than the correlation en- 
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ergy because the kinetic one is a single-particle observable. To this end we can 
use the kinetic equation, of course, an equation which leads to the total en- 
ergy conservation (2). It is immediately obvious that the ordinary Boltzmann 
equation cannot be appropriate for this purpose because the kinetic energy 
is an invariant of its collision integral and thus constant in time. We have to 
consider non-Markovian kinetic equations of Levinson type [1] 



f /.<«> 4 E / j * exp cos [\ (t - f)A,} 

x {fZft(l-fa)(l-K) - /.A(l -£)(!-$}. (3) 
where A E = ^r— + — i2±2l_ denotes the energy difference between 

a 2m a 2m b 2m a 2m b _ 

initial and final states. The retardation of distributions, f a (k,t), f' a {k — q,t) 
etc., is balanced by the lifetime r. The total energy conservation (2) for Levin- 
son's equation has been proved in [6]. 

The full solution of Levinson's equation on the long time scale is a hard prob- 
lem, however, its solution in the short-time region t <C r can be written down 
analytically. In this time domain we can neglect the time evolution of distri- 
butions, fait) = / a (0), and the life-time factor, exp j — = 1. The eq. (3) 
can be integrated with respect to the time and the internal time integral can 
be done. The resulting equation for f(t) represents the deviation of Wigner's 
distribution from its initial value, f a (t) = / a (0) + 5f a (t), and reads 



c p (A sr fdpdq 2 l-cos{^A E } 

X {fati(l-fa)(l-f b ) ~ fah(l-f a )(l-n)} • (4) 

This formula shows how the two-particle and the single-particle concept of the 
transient behavior meet in the kinetic equation. The right hand side describes 
how two particles correlate their motion to avoid the strong interaction regions. 
Since the process is very fast, the on-shell contribution to 5f a , proportional to 
t/r, can be neglected in the assumed time domain and the 5f has the pure off- 
shell character as can be seen from the off-shell factor A E 2 (l — cos {^tA^j). 
The off-shell character of mutual two-particle correlation is thus reflected in 
the single particle Wigner's distribution. 

The very fast formation of the off-shell contribution to Wigner's distribu- 
tion has been found in numerical treatments of Green's functions [4,5]. Once 
formed, the off-shell contributions change in time with the characteristic time 
r, i.e., following the relaxation (on-shell) processes in the system. Accordingly, 
the formation of the off-shell contribution signals that the system has reached 
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the state the evolution of which can be described by the Boltzmann equation, 
i.e., the transient period has been accomplished. 

From Wigner's distribution one can readily evaluate the increase of kinetic 
energy, 

After substitution for 5f a from (4) we symmetrize in k and p and anti- 
symmetrize in the initial and final states which yields the correlation energy 
(2) as 

(«) =-£ / I^M ^y^W - - A). (6) 

This expression holds for general distributions f a . 

Of course, starting with a sudden switching approximation we have Coulomb 
interaction and during the first transient time period the screening is formed. 
This can be described by the non-Markovian Lenard - Balescu equation [7] 
instead of the static screened equation (3). With the same discussion as above 
we end up instead of (6) with the dynamical expression of the correlation 
energy 



ab J ( 2nn > \e(q, 



dkdpdq V£(q) f (k - qf 



2mt, 2mj ' * x 

1 - cos(kA B ) 

x — # — l r a m - /»)(! - h). 



2m a 2m„ 



(7) 



One sees that the bare Coulomb interaction Vc is renormalized by the dielec- 
tric function 



^ M = 1 " ? Vc{q) I (2W %-fi» + iV • (8) 



All internal time integrals are bound to the time dependence of f(t) and in 
the spirit of the above discussion can be carried out. 

To demonstrate its results and limitations, we discuss (6) and (7) for special 
cases that allow for analytical treatment. To this end we use equilibrium ini- 
tial distributions. As the first test, let us evaluate the correlation energy in 



4 



equilibrium which is approached for large times t — > oo. The off-shell factor 
— ^ (l — cos l^tA^j^ then turns into the principle value The equilibrium 
distributions f a are natural for this case. 

At the high temperature limit, where the distributions are non-degenerate 



one can evaluate (6) and (7) for the mixture of particles. We assume the 
plasma consisting of two different types of particles a, b with different masses 
m a ,mf,. Performing a series of integrals we obtain the correlation energies 



ab 



EST M ^ E ((^p) ^ [1 " exp(^)e rf c(6) 

(10) 

The parameter b 2 = (fin) 2 controls quantum corrections. Formula (10) 
is the correlation energy in the second Born approximation of statically screened 
Debye potential as well as dynamically screened one. The latter corresponds 
to known Montroll result in plasma physics [8,9]. 

For identical particles in the classical limit b — > 0, we see that the static 
approximation (6) or (10) underestimates the known value of the correlation 
energy [8,9] while the dynamical result (7) or (10) agrees with this Ward result 



E^t c (oo) = ~(1 - V5F_erfc(6)) = —/uk + o(b), 

(oo) = - e 62 erfc6) = -\e 2 nn + o(b). (11) 

One can see that in the classical limit the static result is just one half of the 
correct Debye-Huckel one (1). The dynamical result yields the correct correla- 
tion energy. This difference can be understood in analogy to the field energy 
of a dipole in an external electric field. If the dipole is already present but has 
just to be ordered, we obtain half of the correlation energy we would have if the 
dipole is formed itself by the field. In our case the static result assumes that we 
have a Debye screening from the beginning before the interaction is switched 
on. The dynamical result counts properly for the fact that the screening has 
to be formed itself which results into twice the correlation energy. 
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In order to compare the time dependency of the correlation energy from 
(6) with molecular dynamical simulations [10], we assume a one component 
plasma which possesses a Maxwellian velocity distribution (9) during this for- 
mation time. From (6) and (7) we find 



dpstatic 
-^corr 



(t) 



e 2 nT , 



dt n 2h 
d EX am (t) e 2 K,T 



dt 



n 



h 



Im 
Im 



(l + 2^)e 2 (l-erf(z)) 



e z i(l-erf(^)) 



2z 



(12) 



where we used z = u p Jt 2 — itj, and Z\ = uo p J2t 2 — it^. This is the analytical 



quantum result of the time derivative of the formation of correlation for stat- 
ically as well as dynamically screened potentials. For the classical limit h — > 
it is easy to integrate expression (12) with respect to times and arrive at 



T7 



static 



(t) = -\e 2 n*\l + ^ - (l + 2u 2 f) exp (u 2 p t 2 ) [1 - erf(^)]}, 



EZT m (t) = -^e 2 nK 1-exp ^ 




1 - erf(^t) 
V2 



(13) 



In Figs. 1 and 2, this formulae are compared with molecular dynamic sim- 
ulations [10] for different values of the plasma parameter T. The parameter 
T = where a e = (j^) 1 ^ 3 is the inter-particle distance or Wigner-Seitz ra- 
dius, measures the strength of the Coulomb coupling. Ideal plasma are found 
for T 1. In this region the static formula (13) well follows the major trend 
of the numerical result, see Fig. 1. The agreement is in fact surprising, because 
we saw that the static result underestimates the dynamical long time result 
of Debye- Huckel (1) *f = v / 3/2r 3 / 2 by a factor of two. 

The explanation for this fact is that we can prepare the initial configuration 
within our kinetic theory such that sudden switching of interaction is fulfilled. 
However, in the simulation experiment we have initial correlations which are 
due to the setup within quasiperiodic boundary condition and Ewald sum- 
mations. This obviously results into an effective statically screened Debye 
potential, or at least the simulation results allow for this interpretation. 

If we go to higher densities of 7 = 0.2 in figure 1 (right) we see that the Debye 
Huckel result is far off the correct equilibrium correlation energy. Nevertheless 
the time scale is still appropriate described by the dynamical result. At still 
higher densities, like T = 0.5 and T = 1, see Fig. 2, non-ideal effects become 
important and the formation time is underestimated within (13). Of course, for 
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Fig. 1. The formation of correlation energy due to molecular dynamic simulations 
[10] together with the analytical result of (13) for a plasma parameter V = 0.1 (left) 
and T = 0.2 (right). The upper curve is the static and the lower the dynamical cal- 
culation. The latter one approaches the Debye-Hiickel result. The exact equilibrium 
correlation energy of MC simulations [11] are indicated by the arrow. 
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Fig. 2. The formation of correlation energy due to molecular dynamic simulations 
[10] together with the result of (13) for a plasma parameter T = 0.5 (left) and T = 1 
(right). The upper curve is the static and the lower the dynamical calculation. The 
latter one approaches the Debye - Hiickel result. The simulation approaches the 
exact correlation energy [11] with higher densities indicated by the arrow. 



higher plasma parameter the Born approximation fails. Nevertheless, formula 
(13) can still almost reproduce the formation time but slightly shorter than 
compared with the simulation. This is due to non-ideality which was found to 
be an expression of memory effects [12] and leads to a later relaxation. For 
strongly coupled plasma the Born approximation fails, of course, to reproduce 
the correct equilibrium value, but reproduces the formation time fairly good. 
The equilibrium value, of course, differs and takes smaller values than the 
Debye - Hiickel result [11,13]. This regime is clearly out of scope of our theory. 



7 



The characteristic time of formation of correlations at high temperature limit 
is given by the time where (13) shows a saturation. This is reached at about 
a time of inverse plasma frequency 



LJ r . 



VthK 



(14) 



The inverse plasma frequency indicates that the dominant role is played by 
long range fluctuation. On the other hand, we also see that the correlation 
time is found to be given by the time a particle needs to travel through the 
range of the potential with a thermal velocity i>th- This confirms the numerical 
finding of [3] that the correlation or memory time is proportional to the range 
of interaction, i.e. rather of short range character. 

In the low temperature region, i.e., in a highly degenerate system /i^$> T, one 
finds a different picture. From (6) follows [14] 



KZ(t) - = e, 



-i low 

v corr 



1 - ^ sin(x) + (|^) 2 (| + [1 sin(x) 



1 _ 1— 1(M.)2 



(15) 



with x = and the equilibrium correlation energy 



l h ' 

arctan — H k 

bi l + bl 



(16) 



We like to point out that the formula (15) 



where the abbreviation is bi = 
is generally valid for any binary interaction and applicable in dense metal as 
well as nuclear situations. The only difference lies in the actual value of the 
equilibrium correlation energy which is dependent on the potantial, of course. 
Indeed, (15) follows from the standard procedure of separating angular from 
energy integrals for low temperatures. The time dependence is carried exclu- 
sively by the energy integrals while the angular averaged interaction results 
in its factor. 



Unlike in the classical case, the equilibrium limit of the degenerate case (15) 
is not reached monotonously but with oscillation that are damped with power 
law t~ l in time. In other words, the correlation energy is rapidly built up and 
then oscillates around the equilibrium value (16). We can define the build up 
time r c as the time where the correlation energy reaches its first maximum, 
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t c = 1.05- 



(17) 



with [i the Fermi energy. Note that r c is in agreement with the quasiparticle 
formation time known as Landau's criterion. Indeed, as argued above, the 
quasiparticle formation and the build up of correlations are two alternative 
views of the same phenomenon. 

The formation of binary correlations is very fast on the time scale of dissipa- 
tive process. With respect to dissipative regimes, the binary correlations can 
be treated as instant functionals of the single-particle distribution and thus 
included into the Boltzmann equation via various renormalizations of its in- 
gredients, would it be the screened Coulomb potential in the scattering rate 
or the quasiparticle corrections. Under extremely fast external perturbations, 
like the massive femto second laser pulses, the dynamics of binary correlations 
will hopefully become experimentally accessible. Even if related measurement 
will not reveal any unexpected features, the experimental justification of basic 
concepts of the non-equilibrium many-body physics is very desirable. The the- 
oretical support to such experiments is mostly based on the non-equilibrium 
Green's functions or the molecular dynamics which both demand expensive 
numerical treatments. Unlike these two approaches, the presented theory fails 
for special systems where the characteristic time of formation of correlations 
becomes longer or comparable with other time scales. For normal systems, it 
provides a simple tool for theoretical predictions. 
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